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Using numerical simulations we study the pinning and dynamics of interacting colloids on periodic 
one-dimensional substrates. As a function of colloid density, temperature, and substrate strength, 
we find a variety of pinned and dynamic states including pinned smectic, pinned buckled, two-phase 
flow, and moving partially ordered structures. We show that for increasing colloid density, peaks in 
the depinning threshold occur at commensurate states. The scaling of the pinning threshold versus 
substrate strength changes when the colloids undergo a transition from one-dimensional chains to a 
buckled configuration. 

PACS numbers: 82.70.Dd 
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Assemblies of interacting colloidal particles in two di- 
mensions (2D) have attracted considerable attention as 
an ideal model system in which various types of equi- 
librium phases and melting transitions can be studied 
conveniently [1] . Some attractive features of this system 
include the fact that the colloid-colloid interactions and 
density can be changed easily and that the individual col- 
loid positions and motions are directly accessible. Fur- 
ther, when a ID or 2D periodic substrate is added to the 
system, unique ordering and melting transitions appear. 
In experiments on colloids interacting with ID substrates 
created using interfering laser beams, a novel laser in- 
duced freezing was observed in which the colloids freeze 
into a crystal as the substrate strength is increased [2]. 
This substrate-induced freezing effect has also been stud- 
ied theoretically [3] and numerically [4]. At even higher 
substrate strengths, a reentrant laser induced melting 
can occur when the laser or substrate strength is large 
enough that the colloids behave one-dimensionally and 
the fluctuations are effectively enhanced, as predicted ini- 
tially in theoretical studies [4]. In this case, the colloids 
form a smectic state. Laser induced freezing and melting 
have subsequently been studied both experimentally [5] 
and theoretically [6,8]. A rich variety of other equilib- 
rium phases for varied colloid density on ID substrates 
have been predicted [8] and in some cases experimen- 
tally confirmed [9]. More recently the crystalline phases 
and melting of colloids interacting with 2D periodic sub- 
strates have been studied and a variety of novel crys- 
talline orderings and melting phenomena were shown to 
occur [10-13]. 

Far less is known about the dynamical interactions 
of colloids with periodic substrates. Recently it was 
shown that dynamical locking effects can occur for col- 
loids driven over 2D periodic substrates when the col- 
loids preferentially move along the symmetry directions 
of the substrate [14-16]. These locking effects may be 
potentially useful as a technique for separating particle 



mixtures. Colloids driven over periodic substrates are 
also a useful model system for studying depinning phe- 
nomena. Similar systems in which depinning on periodic 
substrates occurs include vortices in superconductors in- 
teracting with ID and 2D periodic pinning arrays [17] as 
well as models for atomic friction [18]. 

In this work we consider colloids driven over periodic 
ID substrates, which to our knowledge has not been stud- 
ied previously. We consider parameters relevant to recent 
experiments on colloids interacting with periodic poten- 
tials. For fixed substrate strength and lattice constant, 
we find that as the density of the colloids increases, peaks 
in the depinning threshold occur at various commensu- 
rate fillings. In general we observe a pinned regime, a dis- 
ordered or plastic flow regime, and partially ordered mov- 
ing regimes as a function of driving force. Several types 
of pinned states appear, include pinned commensurate 
lattices, pinned smectics, and buckled pinned configura- 
tions. For high colloidal densities where the colloids form 
a smectic state, we show that as the substrate strength 
decreases there is a change in the scaling of the depinning 
force versus substrate strength when the colloids undergo 
a transition from ID chains to a buckled state. 

We simulate a two-dimensional system of iV c col- 
loids with periodic boundary conditions in the x and 
j/-directions. The overdamped equation of motion for 
colloid i is 
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Here the colloid-colloid interaction force is f* c = 



— SjJi ^iV( r ij)i where we use a Yukawa or screened 
Coulomb potential given by V(rij) — (gjgj/|rj — 
rj|)exp(— K|r-j — r,-|). q^j) = 1 is the colloid charge, 
X/k is the screening length, and ^i(j) is the position 
of particle The substrate is periodic in ID as 

shown in Fig. 1(a), with period d in the y-direction, 
f s = f p sin (2iry/d)y. This is the form expected from the 
modulated laser fields used in the experiments [5]. 
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FIG. 1. (a) Heavy lines indicate the locations of the sub- 
strate minima, (b-d) Colloid configurations (black dots) for 
a system with fixed d, fixed f v = 2.0, and n c — 1.35. (b) 
fd — 0.0 where a pinned smectic state occurs, (c) fd/ fc = 1.1, 
where f c is the depinning threshold. The colloids are moving 
in the y direction, (d) fd/ fc = 1.2. Here the colloids show 
significant triangular ordering. 

The external driving force is fd = fay which could come 
from an applied electric field. The thermal force fr is 
a randomly fluctuating force from Langevin kicks. We 
measure lengths in units of d, the substrate periodicity. 
We have considered various system sizes; here we focus 
on fixed L y = 16d. A triangular colloidal lattice is com- 
mensurate with the substrate for d = V3a/2, where a 
is the colloidal lattice constant. The density of colloids, 
n c , is normalized such that at the commensurate den- 
sity, n c = 1.0. Temperature is reported in terms of the 
melting transition temperature T m for a colloidal lattice 
at the commensurate density in the absence of the sub- 
strate. For the results in this work we focus on the case 
T/T m = 0.5 which is low enough to avoid appreciable 
creep. The initial colloidal state is obtained at zero exter- 
nal drive by cooling the system from a high temperature 
molten state to a lower temperature. A similar procedure 
was used to find the equilibrium states for colloids on 2D 
substrates [10]. The applied drive is then increased from 
fd = 0.0 by small increments. At each increment, we 
measure the average colloidal velocity V y = J^iti v « ' V 
after reaching the steady flow state. 

We first consider the case of fixed substrate lattice con- 
stant d and strength f p and examine the depinning as a 
function of colloidal density n c . In general we find three 
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FIG. 2. (a) Average colloidal velocity V y vs fd/ fp for the 
system in Fig. 1 with f p = 2 and n c = 1.35. Inset: V y vs 
fd/ f P for n c — 0.375. (b) The density of six-fold coordinated 
colloids Pq vs fd/ f v for the system in Fig. 1. 

regimes: pinned, disordered flow, and a partially ordered 
high drive flow. The exact type of order in the pinned 
regime depends on the colloid density. For high densi- 
ties n c > 1.0, the colloids form a pinned smectic state, 
as seen in Fig. 1(b) for n c = 1.35. At most densities, 
including this one, the colloid lattice is not commensu- 
rate with the substrate, so the potential minima capture 
differing numbers of particles and there are dislocations 
present in the lattice with Burgers vectors aligned par- 
allel to the substrate. These dislocations depin before 
the particles do as the drive is increased so the colloids 
generally initially depin into a plastic flow regime or two 
phase flow regime where only a portion of the colloids 
move. In Fig. 1(c) we show the colloidal positions just 
above depinning for fd/ fc = 1-1, where f c is the critical 
depinning force, for the same system in Fig. 1(b). Here 
the colloidal positions are much more disordered, and a 
portion of the colloids are pinned while the other por- 
tion is moving. As fd is further increased, all the colloids 
depin and form a moving partially ordered state where 
the colloids regain a considerable amount of triangular 
ordering, as seen in Fig. 1(d). The colloid lattice also 
realigns in the direction of the drive. 

The different flow regimes can be identified by the 
amount of disorder in the lattice and the characteristics 
of the colloid velocity vs external force curves. In 
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FIG. 3. Flow regions I (pinned), II (plastic), and III (par- 
tially ordered) for fd/f P vs colloid density n c at constant 
f p = 2.0. Open squares: depinning threshold. Filled circles: 
reordering crossover. 

Fig. 2(a) we plot the average colloid velocity V y and in 
Fig. 2(b) we show the corresponding fraction of 6- fold co- 
ordinated particles P & vs the external drive fd/f P for the 
system with n c = 1.35. For a triangular lattice, P@ = 1.0. 
In Fig. 2(a), there are three regions of the velocity force 
curve labeled I, II and III. Region I is a pinned state 
for fd/fp < 0.75. The plastic flow region II occurs for 
0.75 < fd/fp < 0.95 and begins when V y jumps above 
zero. In the plastic flow regime, V y is smaller than it 
would be if all the colloids were moving, and lies below 
the value that would be obtained by a linear extrapo- 
lation of the velocity- force slope at higher drives. For 
fd/fp > 0.94, there is a final jump in V y indicating that 
more colloids are now moving, and the system enters re- 
gion III where V y increases linearly with fd and all the 
colloids move with a uniform velocity. The three regions 
can also be distinguished in P e . In the pinned region I 
where colloids are smectically ordered, there are a sig- 
nificant number of dislocations giving Pq = 0.65. In re- 
gion II flow, a number of the pinned dislocations annihi- 
late; thus, the average coordination number increases to 
Pq = 0.78. In region III, the colloids form a mostly trian- 
gular lattice, giving P 6 = 0.95. In the inset of Fig. 2(a) we 
show the velocity force characteristics for a system with 
n c = 0.375 where the lattice makes a transition directly 
from a pinned triangular lattice to a moving triangular 
lattice. Here the intermediate jumplike features in V y are 
missing due to the absence of a plastic flow regime. 

In Fig. 3 we map regions I through III as a function 
of fd/fp vs n c for fixed f p = 2.0 and fixed d. The de- 
pinning threshold shows several peaks and a prominent 
broad maximum between 0.75 < n c < 1.0, corresponding 
to a triangular colloidal lattice. The peaks occur at 
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FIG. 4. (a) Squares: critical depinning force f c vs f p for 
fixed n c = 1.35. Solid line: f c = f p curve, (b) Buckled state 
at f p = 1.0 and zero drive, (c) Buckled state at f p — 0.75 
and zero drive. 

commensurate densities n c = 0.38,0.26 and 0.185, cor- 
responding to a colloidal lattice constant of a = V2d, 
\/3<i,and 2d, respectively. At these densities there are no 
dislocations present in the system and the colloids pass 
directly from a pinned triangular solid to a moving tri- 
angular lattice as a function of driving force without an 
intermediate disordered flow regime. Region II persists in 
a small sliver between regions I and III near the commen- 
surability at n c = 1 due to the fact that the system was 
not perfectly commensurate and there were a small num- 
ber of dislocations present. There is a broad minimum 
in the depinning threshold centered around n c = 0.45. 
Over this range of density there are significant numbers 
of dislocations present in the system. For n c > 1.0 the 
depinning threshold drops dramatically with density as 
dislocations proliferate and the system enters a pinned 
smectic state. For all densities, once the drive is large 
enough for all the colloids to depin, the system enters the 
partially ordered flow regime III. The region III boundary 
roughly coincides with fd = f p , although for n c > 1.0, 
the onset of region III shifts slightly down in drive for 
increasing n c due to the enhanced colloidal interactions 
at these higher densities. 

The appearance of the three dynamic regimes is similar 
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to what is seen for vortices driven over 2D random disor- 
der [19]. In the vortex case, the plastic flow regime occurs 
when some vortices are trapped in pinning sites while ad- 
ditional interstitial vortices move between the pins. For 
strong driving, the vortices form a moving smectic state 
aligned with the drive, due to the effective transverse 
pinning which is present even at high drives. In the ID 
periodic pinning case considered here, the partially or- 
dered regime at high drives is a polycrystalline moving 
solid rather than a smectic solid since there is no trans- 
verse pinning. 

In experiments with ID periodic arrays, it is very 
straightforward to control the substrate strength by vary- 
ing the laser intensity. Thus we consider the effects of al- 
tering the substrate force f p , and find that when the sub- 
strate is weakened for n c > 1.0, a crossover occurs from 
the ID pinned smectic state illustrated in Fig. 1(a) to a 
buckled configuration. In Fig. 4(a) we plot the critical de- 
pinning force f c vs f p for a system with fixed n c = 1.35. 
The straight line in Fig. 4(a) indicates f c = f p . For 
f p > 1.9 the pinned state is a one-dimensional smectic 
structure. As f p is further increased, f c increases linearly 
with f p and the same pinned smectic state forms. For 
f P < 1-9, patches of the pinned smectic undergo a buck- 
ling transition where the colloids along a single row splay 
out in a staggered manner, as illustrated in Fig. 4(b) for 
f p = 1.0. Here portions of the lattice remain in the ID 
state and coexist with the buckled state. In general the 
buckled portions do not form adjacent to one another due 
to the increased inter-row colloidal repulsion that occurs 
when the buckled state forms. The buckling appears at 
locations of enhanced stress, where dislocations in the 
pattern would be present in the ID pinned case. As f p 
is further decreased the buckled areas grow, as seen in 
Fig. 4(c) for f p = 0.75. The onset of the buckled state 
coincides with a change from the linear behavior of f c 
with f p to the much more rapid nonlinear decrease for 
f p < 1.9. The buckled structure is much more weakly 
pinned than the ID state since the force acting to push a 
colloid over the substrate barrier has an additional con- 
tribution from the staggered colloids on cither side. For 
pinning forces at which a buckled pinned state forms, 
there is both a plastic flow regime and a crossover to 
region III flow as fd approaches f p . 

In conclusion, we have numerically studied the dynam- 
ics and pinning of colloids driven over one-dimensional 
periodic substrates. We find three general regimes, 
pinned, disordered, and partially ordered, which can be 
characterized by the amount of disorder in the colloidal 
lattice and by features in the velocity force curves. We 
map these regimes as a function of colloid density and 
show that peaks in the depinning threshold occur at com- 
mensurate densities where triangular or mostly triangu- 
lar colloid lattices form. As the colloidal density increases 
for strong substrate strengths, the colloids form a pinned 
smectic state which transforms to a buckled structure as 



the substrate strength is reduced. In the pinned smec- 
tic state, the depinning threshold decreases linearly with 
decreasing substrate strength, while in the buckled state, 
the depinning threshold decreases much faster than lin- 
early with substrate strength. 
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